Long-lasting Exponential Spreading in Periodically Driven Quantum Systems 
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Using a dynamical model relevant to cold-atom experiments, we show that long-lasting exponential 
spreading of wave packets in momentum space is possible. Numerical results are explained via a 
pseudo-classical map, both qualitatively and quantitatively. Possible applications of our findings 
are also briefly discussed. 

PACS numbers: 05.45.Mt, 05.45.-a, 03.75.-b, 05.60.Gg 



Apart from idealized exceptions like inverted oscilla- 
tors (I], exponential growth of observable quantities is 
thought to be a rare occurrence in small quantum sys- 
tems. Exponential instability, which is the landmark of 
classical chaos, does not give rise to long lasting, physi- 
cally observable exponential growth, not even in classical 
mechanics, because exponentially fast motion along un- 
stable manifolds is concealed as such manifolds fold in 
homoclinic or heteroclinic tangles. Moreover, it does not 
survive in quantum mechanics beyond a short time scale 
t% — ln(7i)/A 0, where A is the maximal Lyapunov 
exponent. Using a simple dynamical model relevant to 
today's cold-atom experiments, we show in this work that 
exponential quantum spreading (EQS) of wave packets in 
momentum space does exist for a significantly long time 
scale. This is achieved by exploiting two intriguing mech- 
anisms at once. First, our quantum system is set close to 
a quantum resonance condition, in a "pseudo-classical" 
(yet strictly quantum) regime where the quantum dy- 
namics follows trajectories of a fictitious classical system, 
which emerges at small values of a "pseudo-Planck" con- 
stant h* . This "constant" h* measures the detuning from 
exact quantum resonance and can easily be made quite 
small by varying physical parameters, while keeping the 
system in a deeply quantum regime. Second, appropri- 
ate driving fields are chosen such that a pseudo-classical 
stable manifold of an unstable fixed point can support 
easy-to-prepare quantum states. As shown below, these 
factors concur in eventually generating EQS in momen- 
tum space (or exponential energy increase) as a transient 
but long-lasting dynamical phenomenon, whose duration 
is given by (A/h*)\n(C/H*) (A and C are constants for 
small h*) and therefore can be enlarged to a very long 
time scale as parameter H* decreases. The exponential 
rate of EQS scales linearly with H* and is typically small 
for close-to-resonance situations, yet the EQS momentum 
scale can be huge because it scales with 1/h*. Such long- 
lasting EQS does not rest on the chaotic motion which 



is present in the classical limit, but on quasi-integrable 
motion which is instead present in the pseudo- classical 
limit. 

The dynamical model used here is a simple modifica- 
tion of the well-known kicked- rotor model [2| , a paradigm 
of quantum chaos that has been extensively studied both 
theoretically and experimentally. For all previous vari- 
ants of such type of driven systems, quantum spreading 
is known to range from ballistic propagation to super- 
diffusion as well as standard linear diffusion followed by 
dynamical localization Q. Our finding of long-lasting 
EQS is thus an intriguing result in quantum chaos. In 
a broader context, the EQS reported here is different 
from other forms of fast wave-packet spreading like, e. 
the "super-ballistic" diffusion found in a lattice model 
(where the wavepacket spreading can be a cubic function 
of time) and the transient "hyperdiffusion" found in gen- 
eralized Brownian motion [J] (where a square variance 
may increase at the fifth power of time). Furthermore, 
analogous to the novel exponential acceleration of classi- 
cal particles [U , the possibility of long-lasting EQS might 
be useful for designing some robust acceleration methods 
for quantum systems. 

Consider then the following Floquet propagator for a 
double- kicked-rotor model 



U = e-^ T ~ 
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-»TT cos (q) p - iT a 
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with all quantities in dimensionless units. Throughout 
q and p are the canonical coordinate and momentum 
operators. This system was first experimentally stud- 
ied in Ref. |6j by subjecting cold atoms to two period- T 
<5-kicking sequences of an optical lattice potential, time- 
shifted by To < T with respect to each other. Under the 
so-called quantum main resonance condition Th = 4tt, 
achievable in a number of cold-atom experiments 0-[ll| , 
we obtain an "on-resonance double-kicked rotor" model 
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FIG. 1: Two examples of long-lasting exponential quantum 
spreading (solid lines) in RDKR [see Eq. JS)]. (a) K = 5, T = 
1, and Toh = 2tv + h* with h* — 10 -3 . Dotted line represents 
the pseudo-classical result and dash-dotted line represents ~ 
e 2A +* with A+ found analytically, (b) K = 9, T = 1, and 
Toh = 307r/13 + 1.7 x 10~ 4 . Dash-dotted line represents the 
best exponential fitting (shifted) that suggests an exponential 
rate of X' + — 4.1 x 10 -4 . Here and in all other figures, plotted 
quantities are in dimensionless units. 

(RDKR) [l2|,[l3[, described by the Floquet propagator: 

(j R = e -K^-To)^ e -if cos(g) e -iT |^ e -if cos(,j)_ ( 2 ) 

Due to spatial periodicity of the kicking potential cos(g), 
changes in momentum p must occur in integer multiples 
of H. The momentum eigenvalues of the system are hence 
given by mh + /3H, where j3 6 [0, 1) is a quasi-momentum 
variable determined by the initial state and conserved 
by the dynamics. In cold-atom experimental realizations 
of the kicked-rotor model, the initial state is a mixture 
of states with a certain spread of j3, where the spread 
can be made small by loading a Bose-Einstein conden- 
sate into the lattice such that the initial-state coherence 
spans over many lattice constants For j3 = 0, 

the Floquet spectrum of U R plotted against T H € [0, iir] 
yields the famous Hofstadter's butterfly structure [13| . 
with the spectra identical with that of the kicked-Harper 
model [lH for arbitrary irrational values of T h/n [Ta |. 

Representative numerical EQS results for the quantum 
map [/r are depicted in Fig. 1, with the time dependence 
of the mean value (n 2 ) = (p 2 /h 2 ) plotted in a log scale 
for a zero- momentum initial state (hence j3 = 0). As seen 
clearly from Fig. 1, apart from a short transient initial 
stage, log((n 2 )) displays a striking linear dependence on 
the number of iterations of tl R (denoted by t). That is, 
Fig. 1 demonstrates exponential spreading in momentum 
space, with the exponential behavior sustaining over a 
wide range of time and momentum scales. For exam- 
ple, in Fig. 1(a), the exponential spreading occurs from 
t » 10 3 up to t « 10 4 , (n 2 ) increases from HT 1 to 10 6 . 
For the case in Fig. 1(b), the exponential spreading also 
brings (n 2 ) from 10 2 to 10 5 , over a time scale of the order 



of 10 4 . In all cases when such exponential behavior is ob- 
served in our extensive numerical studies, HTq is close to 
2ttM/N, with M and N being odd integers. For instance 
HTo is close to 2ir for the case in Fig. 1(a), whereas for 
the case of Fig. 1(b) it is slightly away from 307r/13. 

To explain the EQS discovered in our numerical 
work, in our theoretical analysis we exploit the above- 
mentioned closc-to-resonance condition. To that end we 
consider cases when HTq is close to a low-order resonance, 
i.e., HT — 2tt + h* with \h*\ -C 1. The reason of the no- 
tation h* will become manifest in Eq. ([3]). We note that 
under the precise condition HTq = 2tt and j3 = we have 
a so-called quantum antiresonance, and that interesting 
aspects of long-time quantum diffusion in the vicinity of 
quantum antiresonance were previously analyzed for a 
different model, describing a kicked charge in a magnetic 
field [l6| . The long-time asymptotic behavior of quantum 
motion is not, however, in the scope of the present work, 
as we are instead interested in the transient, though long- 
lasting EQS. We now proceed by adopting a "pseudo- 
classical" approximation which was successfully imple- 
mented in the study of other close-to-resonance phenom- 
ena observed in cold-atom optics [17]. That is, quantum 
dynamics under close-to-resonance conditions, and hence 
far away from its proper classical limit, may be neverthe- 
less extremely close to the behavior of a fictitious classi- 
cal system, unrelated to the classical limit. Cases with 
(3 = will be studied first. For convenience we also 
set h* > 0. With such specifications, the first factor in 

U n in Eq. ^ reduces to e lT °5t, a time-reversal of the 
normal free-rotor propagator and a feat experimentally 
achieved in Ref. Next we introduce the rescaled 

kicking strength K = Kh* /H and the rescaled momen- 
tum p = ph*/h. As [3 = 0, p is quantized in integer 

multiples of h* , so the relation e ±lT °^ = e ±7iT ( 2 + ^ 
holds true, and using it in Eq. ^ yields 

xe -^(^+^) e -^i?cos(9) _ (3) 

In Eq. ([3]) h* manifestly plays a Planck-constant-like role, 
so that in the "pseudo-classical limit" when h* — >• 0, the 
quantum dynamics described by ([3|) comes closer and 
closer to that of a classical double-kicked rotor whose 
kinetic energy term alternates between ±(^- + 7rp) but 
with a fixed kicking potential Kcos(q). The associated 
pseudo-classical map can then be obtained directly: 

q n+ i = q n + Ksin(p n+ i + q n +i) 
(pn+i + q n +i) = (Pn + q n ) + Ksm(q n ). (4) 

Remarkably, this map can be recognized as a kicked- 
Harper model (l4j . with the canonical pair identified 
as q and p + q. This map has unstable fixed points 
at q = nm, and p + q = nn, with m + n be- 
ing an even number, and Lyapunov exponents X± = 
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FIG. 2: Phase space portrait of the pseudo-classical limit 
of RDKR for K = 5, T = 1, and T h = 2tv + h* with h* = 
3\/2 x KTV P = in (a) and /3 = 4 x 10~ s in (b). 
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In the limit of K — > 0, 



2 ± VK 4 + 4:K 2 j /2 

stable and unstable manifolds of such fixed points con- 
nect in a periodic network, and phase-space is divided 
into cells in the shape of parallelograms, with two sides 
parallel to the q (zero-momentum) axis, and the other 
two inclined by arctan(2). For small nonzero K, split- 
ting of such separatrices occurs, giving birth to an ex- 
ponentially thin stochastic web. Figure 2(a) depicts a 
typical phase space structure, which is indeed a portrait 
deformed from that of a kicked- Harper model [14]. The 
almost straight lines in Fig. 2(a) are the stochastic web. 

The exponential spreading seen in Fig. 1(a) can now 
be understood with our pseudo-classical map. The ini- 
tial state p = corresponds to an ensemble placed al- 
most perfectly astride the stable manifold of the unsta- 
ble fixed point at (q,p) — (tt, 0). The dynamics first pulls 
this ensemble towards this fixed point, and then repels 
the ensemble away from the fixed point at the exponen- 
tial rate A+ along the unstable manifold. A key point is 
that, although this exponential increase involves a max- 
imum change of 2-7T of the pseudo-classical momentum 
p, it translates into exponential behavior of the physi- 
cal momentum p over a huge scale for \h*\ <C 1, because 
p ~ ph/ h* . 

Figure 1(a) shows that results of pseudo-classical and 
quantum calculations nicely agree on the ensemble level. 
It is interesting to note that because the separatrix struc- 
ture in the standard kicked-Harper model [l4| is tilted 
away (45 degrees for vanishing kicking strength) from the 
zero-momentum axis, our exponential spreading mecha- 
nism cannot be generated there, because a state sitting 
on the stable branch of a separatrix would be delocalized 
on a huge momentum scale right at the start. It should 
also be stressed that this mechanism does not rest on 
chaos; although chaotic diffusion inside the thin stochas- 
tic web is itself an intriguing phenomenon [l8j . EQS is 
instead an outcome of quasi-integrability. 

The time scale of the found exponential spreading, de- 




FIG. 3: Numerical t cxp (after which the spreading in mo- 
mentum space deviates from the exponential law) for RDKR 
(black squares) with K — 1, /3 = 0, Toh — 2tt + H* and a vary- 
ing h* . The dashed line represents the theoretical prediction. 



noted by t oxp , can also be estimated. Here we sketch 
a crude heuristic argument, leaving a more formal ver- 
sion for the supplementary material [19|. For a pseudo- 
classical ensemble aligned with the q (zero-momentum) 
axis with a momentum width Sp, exponential momen- 
tum growth p(t) oc <5pexp(A+<) along the unstable man- 
ifold continues until p(t) ~ tt , and therefore t cxp ~ 
X^lniir/Sp). It can be shown that Sp oc h* for the 
zero momentum quantum initial state, so, on account 
of A + K = Kh*/h (for small K), we obtain i cxp ~ 
(h/Kh*) ln(7r/fi*). Figure 3 shows that this theoretical 
scaling with h* is in excellent agreement with the numer- 
ically found time scale of EQS. Additional numerical re- 
sults (not shown) of the X-dependence of t cxp also agree 
with our theory. It is now clear that as h* decreases and 
exact antiresonance is approached, t exp increases and A+ 
decreases [20j |. 

To motivate experimental interest we now examine the 
effects of a non-zero quasi-momentum j3. Let us general- 
ize our previous definition of p, i.e., p = \h*\(p — /3h)/h. 
Still focusing on Toh — 2ir + h* as an example, we find 
the following pseudo-classical map associated with U at 

P'n = Pn + Ksin(q n ), 
q' n = g„+K + ( 2 /3 + lK 

Pn+i = p' n + Ksm(q n ), 

q n+1 = q' n -p n+1 + {2P + l)n. (5) 

Here p' n and q' n are two intermediate variables. Note also 
that this mapping reduces to Eq. ^ if j3 — 0. The phase 
space structure of the mapping in Eq. ([5]) is still peri- 
odic in p with a period 27r, with one example shown in 
Fig. 2(b). Within the regime — tt < p < tt, there are two 
unstable fixed points located at p = 2/3tt, q — 2tt — £ or 
q = tt + £, with £ = arcsin(47r/3/X). As indicated by 
Fig. 2(b), for a sufficiently small (3, the stable direction 



4 




FIG. 4: (Color online) (a) Effects of nonzero /3 on exponential 
quantum spreading. K = 1, To = 1, and Toh = 2ir + h* with 
ft* = 3\/2 x 10~ 4 7r. Solid lines (from top to bottom) are for 
/3/10 -5 = 0, 1.2, 1.5, 1.8 and 5, respectively. Dashed line is 
for /? = f3 CI « 1.688 x 10~ 5 . The initial state is a momentum 
eigenstate with p — j3h. (b) Comparison of the numerical f3 C r 
(black squares) with theoretical /3 CT for h* = 3v2x 10 -4 7r and 
a varying K. 



of the two unstable fixed points is still almost parallel 
to the q axis and the unstable direction is still extended 
in p. The exponential spreading mechanism hence sur- 
vives. The exponential divergence rate can be found from 
a linear stability analysis of the fixed points. A simple 
calculation gives that the rate is now dependent on /3, 



A+(/3) 



In 



K> 2 



K' 2 



K' 4 



/2 



with 



= K 2 - 167r 2 /3 2 . As such, \+(/3) is maximal for 
/3 = 0. This indicates that in the presence of a distribu- 
tion of values of f3 centered at (3 = 0, the spreading dy- 
namics will be dominated by the /3 = component. On 
the other hand, when \/3\ > fi CT = K/4tt = Kh* /{AitK) 
the fixed points disappear, and then the pseudo-classical 
mechanism for EQS is lost. 



In Fig. 4(a) we numerically examine how a nonzero /3 
affects the exponential spreading behavior. There it is 
seen that an exponential spreading law eventually turns 
into a periodic oscillation as /3 increases. Using this tran- 
sition we have also numerically extracted critical /3 val- 
ues as a function of K and then compared them with our 
theoretical analysis of f3 cr . Such a comparison is shown 
in Fig. 4(b), with excellent agreement. The important 
message from Fig. 4 is then the following: to observe 
exponential spreading with nonzero /3 values, it is advan- 
tageous to increase /3 cr , which can be achieved by increas- 
ing the kicking strength K. Assuming that experimental 
uncertainties in /? is around A/3 = 10~ 3 and that the 
near-resonance condition can be achieved at an accuracy 
of h*/(T h) « 1(T 3 , then this requires KT > 4tt, which 
is in the appropriate regime for current cold-atom exper- 
iments dul|. 

As illustrated in Fig. 1(b), EQS is also observed near 
higher-order resonances, i.e., TqH ~ 2ttM/N. There the 
pseudo-classical approach requires nontrivial modifica- 
tions [2lj |. which are not examined in the present work. 
Finally, we briefly mention some interesting applications 
of the EQS found here. First, if we introduce a phase 
shift between the two <5-kicking potentials as in Ref. (22j 
and then tune TqH close to a resonance condition, we find 
that not only (n 2 ) but even the current (i.e. the momen- 
tum expectation value) can increase exponentially over a 
long time scale. A novel type of ratchet accelerator can 
hence be formed. Second, because the dynamical model 
used here can be equally realized in spin chain models 
kicked by an external parabolic magnetic field [12], |23| , 
an exponential propagation of quantum excitation along 
a spin chain is in principle possible. 
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Supplementary Material 

In this supplementary material we elaborate on the ex- 
ponential quantum spreading (EQS) time scale, denned 
as iexp in the main text. The initial ensemble p = is 
almost perfectly aligned with the stable manifold of the 
unstable fixed point at (q,p) — (ir, 0). The ensuing dy- 
namics is produced by the motion along invariant tori 
close to the separatrix [see Fig. 2(a)]: the dynamics ini- 
tially attracts the ensemble to the fixed point and then 
repels the ensemble away from the fixed point at an ex- 
ponential rate. Indeed, these tori close to the separatrix 
can support Floquet eigenstates, and the initial quantum 
state overlaps many of them. 

Our estimate of t cxp consists of three steps. First of 
all, quantized tori are identified with eigenstates of a 
so-called "tilted Harper Hamiltonian" (tHH) (see later). 
This is correct at the lowest perturbative order. Let |0) 
denote the state of (angular) momentum 0, i.e., the ini- 
tial quantum state. It overlaps the eigenstates of tHH 
whose energies lie in a window of (-5E, 0). An estimate 
for 8E is derived by computing the spread of tHH in the 
state |0). As detailed below, 5E is given by: 

5E ee v/<0|i% t |0> - <0|tf tllt |0> 2 = \kh* . (6) 

Second, we estimate the classical half-period T(E) of mo- 
tion along tori in the vicinity of the separatrix. To this 
end we shall use the Kicked-Harper Hamiltonian in con- 
tinuous time. This is correct at the lowest perturbative 
order and the separatrix is at energy E = 0. We then 
obtain T(E) for E < and \E\ < 1. It will be shown 
that T(E) diverges logarithmically as the separatrix is 
approached, and at the leading order 

T(E) ~ R- 1 \tl{~2K / E) . (7) 

Finally we identify i cxp with T(c5E), where c is an ad- 
justable fitting parameter. Hence we arrive at the follow- 



ing scaling discussed in our manuscript: 

t oxp ~ K-Hni^-) ~ -A-ln(A) (8) 
\ft*> Kh * \ h *> w 

with A = 4/c. Next we sketch the calculations leading 
to Eq. © and Q. 

Energy width of initial state - The quantum double 
kicked rotor map U^ =0) = V^Uo V U with 

may also be written as 

where 

H = Kcos(q), 

Hi = V- 1 H 2 V. (10) 

A Floquet Hamiltonian can be defined as a self-adjoint 
Hf such that exp(—iHf/h*) — U^~°\ None is explicitly 
known about Hf. However for small kicking strength an 
approximation is provided (via Trotter-Kato) by 

Hf « H tm = H + Hi, 

which may be termed the "tilted Harper Hamiltonian" 
(tHH). Its pseudo-classical limit is obviously 

H tm = K[~cos(p + q) + cos(g)] , (11) 

however it does not really coincide with the back- 
quantization of it. Its matrix elements in the eigenba- 
sis \n) (n £ 1*) of the angular momentum p are easily 
computed noting that 

^° = \ R E(l n >< n+1 l + \n){n-l\), 

n 

whence, using Eq. (fl"0|) . 

Hi = --K ^{ e - ift *(™+V2)| n )( n + 1 | 

n 

+ e ifi*(n-l/2)| n ^ n _ 1 || j ( 12 ) 

and so 

<n|# ult |0) = ~(l-e^/ 2 ) *(n+l) + \{l-j n 'l 2 )5{n-l). 

Using this result, we have (0|7J tilt |0) = and 
= E|(0|i? tiIt |n)| 2 

= 2 J fC 2 [l-cos(/f/2)] ~K 2 tf 2 /A. (13) 

Period of classical motion near separatrix - The clas- 
sical period of motion along a tori close to the sepera- 
trix was long known to be diverging as the initial condi- 
tion approaches the separatrix For self-completeness 
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we give a detailed derivation here. For convenience we 
consider the following Hamiltonian, which is canonically 
equivalent to the classical tHH of Eq. ([TT]) : 

H = K [ cos(q) — cos(p) ] . 

The action variable at H = E with E < is 

1(E) = ^fpdq 
1 f n 

= - dq arccos[cos(g) - EK^ 1 ] (14) 

Jan 



where cos(go) = 1 + EK 1 and so , as E 

q a ~ V-2EK- 1 . 
Half-period at energy E reads as 

dl(E) 



0, 



T(E) = 7T 



dE 



K~ 



dqo 
dE 

dq 



arccos[— EK 1 + cos(<7o)] + 
1 



1 - {-EK- 1 +cos(g)) 2 



(15) 



The first term vanishes, and the last integral splits as 
follows: 



dq 



TV — go 



dq (...); 



Here the second integral is ~ (ttK)- 1 as E — > 0, hence 
at leading order will be neglected. As to the first inte- 
gral, expand at the first Taylor order in E the integrand 



becomes 



1 - [cos(g) -EK- 1 ]' 



a/I - cos 2 (g) 
E 1 



2K 



^[l-COS^g)] 3 

(16) 



Substituting this into the integral, the second term is 
found to yield a contribution which does not diverge in 
the limit E — > and so will be neglected. The first term 
yields a contribution that rewrites as 



K^ 1 / dq csc(g) 
J qo 

7r/2-(Jo/2 

R- 1 ln[t&n(q)] 



qo/2 



-IK' 1 In 



tan(A/- J B(2 J ft:)- 1 



- K^M-IKjE) , (17) 

and so, putting all the above together, we have (at the 
leading order) the following 



T(E) - K-Hn^K/E) , 



which is exactly the Eq. (7) we need. 



(18) 
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